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1. Introduction: Motivation 

The phase diagram of strongly interacting matter represents the, perhaps, most challenging 
open problem within the Standard Model of elementary interaction. Its most prominent feature, the 
deconfinement transition line between hadrons and partons, has been first addressed by R. Hage- 
dorn [|l[], well before the advent of QCD, in his studies of the limiting temperature occuring in 
hadron-resonance matter. The resulting phase boundary, at about T = 160— llOMeV (that concurs 
with an energy density of about 1 GeV / fwr'), was subsequently understood [^, |3|] as the location of 
the QCD hadron to parton deconfinement transition. At such low temperature, and corresponding 
mean Q^, deconfinement can not be the consequence of QCD asymptotic freedom - the perturbative 
QCD mechanism that was first envisaged iQ, - but falls deeply into the non-perturbative QCD 
sector, as does the resulting confined hadron structure. In fact, non-perturbative QCD theory on the 
lattice has, for two decades, postulated that the hadron-parton phase transformation occurs in the 
vicinity of Hagedom's limiting temperature [§, 0], at zero baryochemical potential. 
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Figure 1: Sketch of the QCD phase diagram in grand canonical variables T, pLg- A critical point E is 
indicated which ends the crossover transition domain at /ig below E. Also shown are the hadronic chemical 
freeze-out points resulting from statistical model analysis, from RHIC down to SIS energies. 



Conditions reached in the cosmological expansion evolution, and closely approached at RHIC and 
LHC energies, in collisions of heavy nuclei. Subsequent, recent developments in QCD Lattice the- 
ory have overcome the technical limitation to the case of zero baryochemical potential [Q, ^ p^ ], 
with extrapolations of the deconfinement phase boundary, upward to about pLs = 500 MeV in the 
{T,}JLb) plane, including hints of a critical point of QCD. This domain coincides with the energy 
region of A + A collision study at CERN SPS energies, 5 < < \1 GeV , which has created several 
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Figure 2: Light quark number density fluctuation (susceptibility) from a two-flavour lattice calculation 
for quark chemical potentials Hq = 0.75 and 150 MeV, respectively. 



non-trivial, much discussed experimental observations to which we shall turn below. The expec- 
tations raised both, by the recent lattice QCD developments, and by the former SPS experiments, 
concerning the physics of the phase boundary at non-zero baryochemical potential, have recently 
been reflected in new experimental efforts, notably with the plans for RHIC running at such low 
energy, and CERN re-running with SPS experiment NA61, and with the FAIR facility developed at 
GSI. The present article will deal with the corresponding physics goals, and challenges. We begin 
with a short recapitulation. 



Our present view of the QCD phase diagram is sketched in Fig. It is dominated by the con- 
finement boundary line above which some kind of parton plasma should exist. Up to about 
500 MeV it results from recent lattice QCD extrapolations [M 



1Q| ]; its further extension 



13]. Lattice studies agree 



stems from "QCD-inspired" theoretical considerations [ ]11[ , [1 
that the domain toward Hb = should feature a rapid crossover transition, with a coincidence of 
deconfinement and QCD chiral symmetry restoration. As the deconfinement transition at high Hb 
is generally expected [[Tl|, 12, 13] to be of first order one is led to conclude on the existence of a 
critical point E, at intermediate /Xg. In fact lattice QCD has reported first evidence for the existence 
of a critical point [|, g [I|], to which we shall turn below. Fig. |l| also shows experimental in- 
formation concerning the {ji.B,T) values at which the hadron-resonance species population freezes 
out from the dynamical evolution, in A + A collisions of mass A « 200 nuclei at various ener- 
gies ranging from top RHIC, via SPS and AGS, down to SIS energy. Remarkably, this freeze-out 
occurs in "hadro-chemical" species population equilibrium, as reflected in the Grand Canonical 
Gibbs ensemble of the statistical hadronization model [16, 17, 18 1. We observe that the freeze-out 
points merge with the lattice QCD phase boundary, from top SPS to top RHIC energies. It thus 
appears that hadron-resonance species freeze-out coincides with the QDC confinement conditions, 
leading to hadronization, at low }JLb- In fact, the idea that the phase transformation "gives birth" to 
an equilibrium hadron-resonance population was first formulated by Hagedom [19] and has been 
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further explored recently [20, 21, 22, 23]. We would, thus, locate experimentally the parton-hadron 
phase boundary at /I5 —>■ 0. However, Fig. [T] shows that the freeze-out points fall well below the 
conjectured confinement boundary as /ig — > 500 Me V. It remains an open, challenging question to 
identify a transport mechanism, active during expansion after confinement, which could deliver the 
system in species equilibrium at much lower temperatures. We turn to this topic in the next section 
but note, for now, that the recent hypothesis [ p4| ] of a further QCD phase (labeled "Quarkyonium"), 
occuring in between confinement and hadronization, might provide for another phase transition, at 
high Hb, which enforces hadro-chemical equilibrium conditions at freeze-out. 



Early consideration of a critical point at finite /Xg arose from quark susceptibilities in lattice QCD 
[^5|], as illustrated in Fig. ^ The calculation employed two dynamical quark flavours and showed a 
smooth ("crossover") behaviour at the critical temperature Tc for Hb = ^ IJ^q = 0, but a dramatic sus- 
ceptibility peak developing toward /Xg = 3 = 450MeV. Such a quark number density fluctuation 
was implied, as a characteristic consequence of a critical point, in the pioneering 1999 article of 
Stephanov, Shuryak and Rajagopal p^]. The most recent lattice results, of a 3 flavour calculation 
with near-realistic quark masses [^] are shown in Fig. ^. The sharp peak structure in the light 
quark sector remains, at /ig = 2 Tc, and even the strange quark susceptibility exhibits a (dampened) 
irregularity. 

Further toward gaining hints for the existence of a critical QCD point, from state-of-the-art lattice 
calculations, C. Schmidt of the RBC-Bielefeld collaboration showed [28] at QM08 an analysis of 
the radius of convergence, observed in the Taylor expansion coefficients (toward high Hb) of the 
pressure, as illustrated in Fig. ^. By definition, the Taylor expansion convergence ends at a critical 
point. In fact, the second order estimate of the radius of convergence exhibits a first indication of 
a nonmonotonic decrease at T « 195 MeV, /ig 500 MeV in this calculation. Higher orders in 
the pressure expansion are clearly required to pin down this first result - with obvious implications 
concerning lattice QCD computing power advances. 

The energy domain, of interest to the search for manifestations of a critical point, and the adjacent 
onset of a first order nature of the phase transition, is broadly outlined by such lattice observations: 
to fall into the interval 250 < Hb < 500 MeV. This corresponds, almost exactly, to the SPS range 
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Figure 4: Radius of convergence of the Taylor expansion coefficient for the pressure, in the T, iib plane, 
from a lattice 3-flavour calculation indicating non monotonic structure at jis ~ 500 MeV. 



6.5 < < 17.3 GeV , also including the top energy reached in former AGS investigations. At 
the low end of this interval, we might expect the "onset of deconfinement" [|^]: the dynamical 
trajectory of central heavy nucleus collisions would reach its turning point (between initial com- 
pression/heating and final expansion/cooling) in the close vicinity of the parton-hadron boundary 
of Fig. [1} 

Moreover, the dynamical time scale is the (relatively) slowest at the turning point, of order several 
fm/c at these energies, such that equilibrium might be approached and the evolution would exhibit 
a sensitivity to thermodynamic features such as a critical point, and to bulk matter phenomena like 
a first order phase transition. Please note that the A + A dynamical evolution never proceeds on 
a "line" in any phase diagram, nor does it collectively reach any "point" in the {T,}JLb) plane, due 
to the relatively unfavourable surface to volume ratio creating density gradients in the "fireball" 
system - the evolution proceeds within a relatively broad band of local conditions. One would thus 
not expect a sharply singular behaviour of any observable with y/s. 

Seen in this hght, the relatively sharp peak observed by NA49 in the 1%^ excitation function 



[30] is remarkable. Fig. || shows the dependence of this ratio, from threshold behaviour (at 
AGS energy) to peak at low SPS energy, with decline toward top SPS and RHIC energy. This 
is probably the sharpest dependence on that could be expected. Is it a signature of "onset of 
deconfinement"? We will turn to this open problem below, and in the next section, but note here 
that Fig. ^ also illustrates the phenomenon called "strangeness enhancement in A + A". Note that 
represents the bulk of (anti-)strangeness production, and that the /n'^ ratio thus approaches 
the global strangeness to entropy ratio. Fig. ^ also illustrates the data for elementary minimum bias 
p + p collisions which exhibit no peak structure. The dramatic overall increase in central A + A 
collisions has been understood as a consequence of the formation of a large, quantum number 



coherent interaction volume in A + A collisions, absent m p + p interactions [31]. Within the 



terminology of the statistical hadronization model, this effect is reflected in a transition from a 



canonical to a grand canonical Gibbs ensemble description [31], where "canonical strangeness 
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Figure 5: Energy dependence of the / ratio in central A+A and in p + p collisions, from AGS to 



RHIC energies [22 , 301 



suppression" in elementary, small volume collisions is seen to arise from strict quantum number 
conservation constraints at the microscopic level, fading away in large, coherent volume solutions, 
corresponding to the A+A case, which lead to strangeness saturation. The / 71+ signal in Fig. || 
thus tells us that A + A collisions create a large quantum number coherent "fireball" volume before 



hadronic freeze-out sets in [32]. 



The K/n ratio should also be "robust" insofar as it will survive final state interaction during the late 
stages of hadronic expansion, during which strangeness and entropy are approximately conserved. 
This motivates the analysis of its fluctuation at the event-by-event level [33, ^]. Within finite 



sub-intervals of rapidity space this might be caused by the lattice quark density fluctuations shown 
in Fig. ^ As shown in Fig. ^(left) the initial measurement [33] of the quantity Odyn{^) gave a very 



small result, the signal being almost fully exhausted by the mixed background. Fig. ^ (right) shows 
that this remains so at RHIC energy [Q] but that Odyn exhibits a significant rise toward low SPS 
energy [35], where — > 450MeV, the region of interest from Figs. ^ and^ A simulation with the 



UrQMD transport model shows no such dependence. Randrup and collaborators have suggested 
[36] to approach, both the ensemble average K/n (= + /n^ + n^) ratio, and its event- 



wise fluctuation, not as a critical point consequence but as effects of the spinodal decomposition 
mechanism that could develop during a first order phase transition. The global grand canonical 
fireball (constrained by charge and strangeness conservation on average) would decompose into 
separately expanding "bubbles" in phase space which could carry nonzero net strangeness. Fig. ^ 
shows that the ensemble average kaon yields (left) exhibit no difference from the grand canonical 
global values (as could have been expected from the percolation models for strangeness one [p^]) 
but that the event by event multiplicity variance does indeed grow well beyond the GC value. 
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Figure 6: Left: Event by event fluctuation of the charged kaon to charged pion multiplicity in central 
collisions of Pb + Pb at -^5 — 17.3 GeV, compared with background from mixed events (histogram) p3[|. 
Right: energy dependence of adyn for kaon to pion fluctuation at SPS and RHIC [ p4[ ] energies, compared 
to UrQMD prediction. 



toward high similar to the data in Fig. ^. 

We shall return to these observables in the next sections but end, at this point, our introductory tour 
de horizon of the field - as is stood at the previous Firenze CPOD workshop We now turn to 
the present physics of the QCD phase diagram. 



2. Recent Developments 

2.1 The critical Point Focusing Effect 

Asakawa and Nonaka have developed a hydrodynamical model which implements a critical re- 



gion [ )38[ , |39p, causing a sensational attractor or focusing effect on the dynamical trajectories of the 
system expansion at different ^/s. Their work is illustrated in Figs. |8| and |^ Fig. || (right) gives 
dynamical trajectories in the absence of a CEP but including a first order phase transition taking 
place at the parton-hadron boundary line. 



It is important to recall that in reality an A + A central collision creates, not such lines but relatively 
broad bands in the [T, Hb] plane but we infer, nevertheless, that the regions in which the equilibrium 
period of the fireball expansion enters the figure (from above) are monotonously shifting (upward 
and) to the left with increasing ^/s. The influence of the confinement mechanism considered leads 
to back-bending and some further reshuffling (characteristic of the phase transformation type im- 
plied) of trajectories. There is no hadronic freeze-out considered in this model. But we know from 
Fig. [T] that the data are arranged in the close vicinity of the coexistence line, falling off at large 
Hb- It is important, now, to note that each expansion trajectory has to pass through such a chemical 
freeze-out "point". So, from this version of the hydro model we infer a smooth, monotonous rela- 
tionship between ^/s and freeze-out points, and this is, in fact, the result of the systematic analysis 
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Figure 7: Effect of hydrochemical spinodal decomposition on kaon production in central A+A collisions 
[Q. Comparison to a single grand canonical fireball is made, with increasing Hb, for the average multiplicity 
(left) and for the multiplicity fluctuation width. 



with the statistical model [[T^, 17, [T^, ^ But we also learn that this relationship must be 
sensitive to the position of the coexistence line, and to the nature of the phase transformation. 
This is true, in particular, if a CEP focusing mechanism is introduced, as illustrated in Fig. |8| (left). 
The trajectories enter, at model-time zero, at T = 200 MeV and Hb = 145, 461 and 560 MeV, 
respectively, similar to the case with no CEP, but get substantially defocused. The relationship 
between y/s and freeze-out /I5 would be dramatically affected. This is further illustrated [ |39| ] in 
Fig. ^ which shows the expansion trajectory in the close vicinity of an entropy to baryon number 
density ratio (conserved during hydro-expansion, thus a "label" of the individual trajectory) of 
about 25. The same chemical freeze-out point is reached by 3 different trajectories corresponding 
to different ^/s and assumed mechanism for crossing of the phase boundary: a first order transition 
without (FO) and with (QCP) presence of a CEP, and a crossover transition. 
Without turning to discuss the degree of realism of this model (please consult [38, 39]) we conclude 
that the freeze-out points in Fig. |I| need to be revisited. 



2.2 Hadrochemical Freeze-out at High Hb 

Turning to the hadrochemical freeze-out points in Fig. |l|, we note, first of all, that the statistical 
Grand Canonical ensemble of hadrons and resonances can not, in its ideality, be referred to the re- 
alistic situation encountered in central A+A collisions without substantial additional assumptions 
and approximations |22|]. Think e.g. of rapidity windows, surface "corona" effects, canonical 
strangeness suppression, finite hadron size, etc.. Thus, each of the various "schools" of statistical 
model analysis has, over time, changed (improved) its technical details, and, moreover, different 
choices of dealing with such second order concerns have been persued. Some scatter of the pub- 
lished [r, n] freeze-out points for each ^/s is the consequence and, taken together, they convey 
the impression of a band rather than a curve; this band has then been well described by a smooth 
parametrization [23, 4C, 41], the "freeze-out curve". 

However, toward our present concern about non-smooth excursions from the freeze-out curve, 
we take the opposite approach. We illustrate in Fig. IC exclusively the freeze-out points from a 
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Figure 8: Left: effect of a critical point on the hydrodynamic expansion trajectories in central A +A colli- 



sions at various energies [38 1. 

Right: the same initial trajectories without a critical point but showing the effect of a first order phase tran- 
sition in the EOS. 



single most recent systematical investigation of all available data [23]. The overall trend remains 
but excursions are obvious, perhaps significant in view e.g. of the CEP trajectory re-focusing 
discussed above. These results are the consequence of an attempt to describe the K^/n^ peak 
of Fig. H that we illustrate in Fig. |ll[ It is the novel ingredient of analysis [p3|], effective here, 
to consider more explicitly the energy dependence of the hadronic resonance spectrum that enters 
the Grand Canonical hadron plus resonance Gibbs ensemble. The model now captures the decline 
in the K^/n^ and A/n^ signals at ^/s above about 8GeV, the sharp peak structure in K^/n^ 
still elusive. However, it is also obvious from Fig. 1 1 that the data in the SPS energy range need 
significant improvement. A far higher accuracy is also of crucial importance to our next topic. 

We return to the hadronic freeze-out points, smooth or not in Fig. |l| and |l^, and focus at their fall-off 
from the confinement line suggested/conjectured by lattice QCD. Note that the term freeze-"out" 
refers to the Gibbs ensemble of on-shell hadrons and resonances in vacuum, as it is employed in the 
statistical model of hadronization. There is no mean field, and zero virtuality. The hadron gas is still 
interacting due to the presence of the full resonance spectrum, but the number distribution of species 
(after resonance decay) stays constant (frozen-in) from this point on. Note that T = 160— llOMeV 
represents a relatively cold medium, with small average inelasticity. This is illustrated in Fig. 12 
(top) by a hydro plus UrQMD calculation by Bass and Dumitru [ p2| | for central Au + Am collisions 
at top RHIC energy. The yield distribution after Cooper-Frye termination of the hydrodynamical 
evolution at T = llOMeV is compared with the outcome if a UrQMD hadron-resonance cascade 
expansion phase is attached. It causes only minor changes of the population, which thus stays 
frozen-in, for the case of hadronization of T = 155 MeV and jJ-s ~ 20 MeV. A similar observation 
is made at SPS energy [43]. Fig. 12 (bottom) illustrates the same procedure at y(s = 8.7 GeV for the 
case of hadronization at l50MeV and /ig = 350MeV. These calculations illustrate the well-known 
fact that the chemical relaxation time of a binary hadron-resonance cascade mechanism is far too 
long [ 44 ] to achieve equilibrium or, turning around, to abolish an equilibrium previously installed 
by the hadronization process [20, pT, 



I 45]. 



We are led to conclude on an apparent impasse: if confinement and hadronization occurs at the 
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Figure 9: Same model as Fig. 8 p^: three initial expansion trajectories characterized by their entropy to 
baryon density ratio S/ns, experiencing the effect of a critical point, a first order, and a cross-over phase 
transition, respectively, end up in a common freeze-out point. 
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phase boundary, also at higher Hb, there is no simple transport mechanism that would carry the 
chemical equilibrium from hadronization at, say, 150 — l60MeV downward to the observed hadronic 
freeze-out equilibrium at, say, 135 — 145 MeV. 

One way out: the data entering the statistical model analysis are insufficient at < 9 GeV, jls ^ 
300 MeV, because the phi-meson and higher hyperon multiplicities are inaccurate, or missing. In 
other words, chemical equilibrium might in fact be incomplete or missing, the freeze-out points 
at high /ifi thus appearing illegitimately on the equilibrium phase diagram of Fig. [l| However, 
there are no indications of such a deviation from equilibrium as far down as y/s = 8.7 GeV, pLs = 
3 80 MeV, as is shown in Fig. [13] by the fit pi] ] to the corresponding SPS, Pb+Pb central collisions 
data of NA49. Lower energies will hopefully be accessible during the planned RHIC runs at low 
energy. 



The alternative: there should be a new state of QCD intervening between the color confinement and 
the hadronic freeze-out curves. This conclusion is also hinted at by the observation [^] that the en- 
ergy density along the lattice-inspired color confinement line far exceeds the canonical 1 GeV / fwr' , 
at }JLb —>■ 500 MeV. Not a likely hadronization condition. So that it is tempting to identify the 
freeze-out points as the Hagedorn limiting temperature p6[ p7| ], i.e. as the limit of existence of 
hadron-resonance matter, at high /ig. Finally, there does actually even exist a candidate of an in- 
tervening, new QCD phase, conjectured by McLerran, Pisarski and collaborators [24], and named 
"quarkyonic matter". Fig. |T^ gives a sketch of the corresponding phase diagram of QCD, which 
features a phase boundary of color confinement and a further one of hadronization from quarkyo- 
nium to hadrons plus resonances (about which nothing is known in detail as of yet). This would 
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Figure 10: The QCD phase diagram of Fig. 1 but with statistical model hadron freeze-out points from an 
analysis pa] including an energy dependence of the hadronicresonance spectrum. 



finally explain the origin of the chemical equilibrium, implied by the hadronic freeze-out points 
at high Hb- to closely coincide with the latter phase boundary ||2^]. In keeping with the observa- 
tion [20, ETl E2|, Bsl] that hadrons do not dynamically acquire chemical equilibrium but are "bom 



into it" by a phase transition which enforces phase space dominance via the emerging 



hadron-resonance spectrum [48]. In this view, the apparent equilibrium in A+A freeze-out re- 
quires synchronization by a global change of degrees of freedom [p8[], unlike freeze-out from a 
smooth dilution process caused by expansion, such as the cosmological synthesis of light nuclei, 
or the final momentum space decoupling of diluting hadronic cascades, which necessarily lead to 



sequential freeze-out [|22p, in order of total cross sections. 

It will be important to measure complete hadron species multiplicity distributions, down to the 
anti-Omega, at energies as low as = 6 GeV (where this multiplicity is as low as 0.02). A 
challenging experiment task, given the low luminosities expected here, at RHIC, but the only test 
of the quarkyonium hypothesis known thus far, as it addresses the domain of Hb —>■ SOOMeV where 
the present, preliminary freeze-out points fall really far below the lattice QCD confinement Une 
estimate (Fig. |l|). 



3. Fluctuation Signals from Event-by-Event Physics Analysis 

Central collisions of heavy nuclei create bulk hadrons, at pr < 2 GeV, with a multiplicity 
per event that is high enough to perform analysis of physics observables at the level of single 
events. Experiments such as NA49 at the SPS, and STAR, PHOBOS at RHIC were designed. 
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Figure 11: Energy dependence of the /n^ and A/;r+ ratio with statistical model description [ ^3[ ] as in 
Fig. 10. 

specifically, for this task. Most copiously produced charged hadrons are and K~ , protons (and 
antiprotons at RHIC), and charged pions, with recorded average multiplicities per event ranging 
from about 5 for at low SPS energy, to several hundreds of pions at the top RHIC energy. Each 
central Pb + Pb or Au + Au collision thus represents, in statistical model terminology, a "grand 
micro-canonical ensemble", where the term "grand" indicates the large production volume, but the 
term "micro" reminds us of the fact that each individual large volume event must overall strictly 
conserve electric, baryonic and strangeness charge conservation. This complicates the seemingly 
straight-forward idea of event by event analysis, as far as conserved quantities are concerned [p^]. 
The first considerations of fluctuations, occurring as a consequence of a QCD critical point [p^], 
referred to non-conserved quantities: eventwise analysis of multiplicity density, and transverse 
momentum fluctuations. We shall turn to the latter observable below, then to continue with hadron 
ratio and higher moment fluctuations, such as kurtosis, which make contact to the recent lattice 
QCD calculations. 

3.1 Transverse Momentum Fluctuations 

The classical, pioneering article of Stephanov, Shuryak and Rajagopal considered a crit- 
ical point of QCD within the framework of the chiral symmetry restoration transition of QCD. 
Complementary to the deconfinement transition, between bound color singlets and finite color mo- 
bility, this transition addresses the transition between massive hadrons (with mass caused by non- 
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Figure 12: Top: effect of the hadronic expansion period on the hadro-chemical composition established at 
hadronization of a hydrodynamic partonic expansion, for central Au +Au at top RHIC energy [p^. 
Botton: same comparison for Pb + Pb at SPS energy 40 AGeV ( = 8.7 GeV). 



perturbative parton condensates in the non-perturbative vacuum) and approximately massless light 
quarks. At low /ig these two QCD transitions appear to occur synchronously [ |l4] , |l5|]. The order 
parameter of the chiral phase transition is the mass of the sigma field meson which approaches zero 



at the critical point, causing "critical" field fluctuation [26]. The linkage to experimentally observ 



able quantities consists in the o ^2n decay near threshold, causing charged particle multiplicity 



density fluctuations [50] as well as transverse momentum pair correlation, and overall transverse 



momentum fluctuation, in the vicinity of a QCD critical point. 

These predictions have inspired numerous experimental investigations, first begun at SPS energy. 



Fig. 15 shows [51 ] the distribution of the eventwise average charged particle pj in central collisions 
of Pb + Pb at top SPS energy, y/s = 17.3 GeV . A perfect Gaussian. The histogram represents the 
signal derived from randomized Monte Carlo sampling of tracks throughout the recorded event 
ensemble: the "mixed background". The real event data show no significant deviation from the 
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Figure 13: Statistical model fit p4| ] to hadronic multiplicities obtained in central Ph + Pb collisions at 

40 GeV ( = 8.7 GeV) by NA49. 



mixed background, over four orders of magnitude: a first "nuU"-result, at this energy. 
We shall see in the next section that such signal comparisons to a mixed background can also lead 
to significant differences, for other observables. But we note, for now, that this mode of analysis 
may be intuitively convincing but leads to considerable difficulty, as far as quantitative compari- 
son to theoretical models is concerned. The mixed background distribution reflects non-trivially 
convolution of experimental biases stemming from acceptance, tracking and particle identifica- 
tion efficiencies, momentum resolution limitations etc., which also might affect the single event 
observation in a way different from the ensemble averaged background, due to impact parameter 
fluctuation, and to event-by-event efficiency variation, correlated e.g. to the eventwise multiplicity 
density distribution. Thus, several alternative measures of fluctuation have been suggested. One 
can base a signal of dynamical pT fluctuation on the binary correlation of particle transverse mo- 



menta in a given event i.e. on the covariance ipriPTj) of particles /, j in one event []52|]. Of course, 
the covariance receives contributions from sources beyond our present concern, i.e. Bose-Einstein 
correlation, flow and jets (the jet activity becomes prominent at high ^/s, and will dominate the pr 
fluctuation signal at the LHC). In covariance analysis, the dynamical pj fluctuation (of whatever 
origin) is recovered via its effect on correlations among the transverse momentum of particles. 



Such correlations can be quantified employing the two-particle pr correlator p3[ ] 

1 



{ApTiApTj) 



M 



n N(k) N{k) 
pairs )t=l ,-=1 j=i+l 



(3.1) 
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Figure 14: Sketch of the QCD phase diagram including the hypothetical quarkyonium phase |24|. The 
hadronic freeze-out points now follow the hadronization line. 



where Mpdrs is the total number of track pairs of the events k contained in the entire ensemble 
of n events, A^(^) is the number of tracks in event k, and \pTi = pn — Pt where pj is the global 
ensemble mean pj. The normalized dynamical fluctuation is then expressed as 



a{pT) 



dyn 



{ApTiApTj)/PT 



(3.2) 



It is zero for uncorrelated particle emission. 



Fig. |16| shows the analysis of pj fluctuations based on the pr correlator, for central Pb+Au SPS 



collisions by CERES [ ]53[ ] and for central Au +Au at four RHIC energies by STAR []52|]. The signal 
is at the 1 % level at all ^/s, with no hint at critical point phenomena. Its small but finite size could 
arise from a multitude of sources, e.g. Bose-Einstein statistics. Coulomb or flow effects, mini-jet- 



formation, but also from experimental conditions such as two-track resolution limits [53]. NA49 



has employed [54| the measure ^(pr) defined as 



^{pj 



<Z2> /= 

\/ z 

<N> ^ 



(3.3) 



where Zi = pn — PT^ with pj the overall inclusive average; for each event one calculates Z = 
with N the multiplicity of the event. With the second term the trivial independent particle emission 
fluctuation is subtracted out, i.e. vanishes if this is all. 



Fig. ^ shows the NA49 result 1 54 ] for central Pb + Pb collisions, and indeed is compatible 
with zero throughout. The SPS energies, from 15S AGeV to 20 AGeV, are represented here by 
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Figure 15: Distribution of event average pj in central Pb + Pb collisions at =17.3 GeV, compared to 
mixed background histogram . 



the prevailing baryochemical potential values [41]. The figure includes a study of the effects that 



might be expected from a critical point, the position and width t, of which are taken from the lattice 
calculations ^]. The expected signal amplitude is then calculated [58] in the framework of 
ref. [|^]. It is relatively small for realistic t, but the data do not allow for a final conclusion; at 
least there is clearly no extremely long range of correlation, such as 6 fm. We see that a substan- 
tial effort is required to conclude on this challenging critical point search proposal. Furthermore, 
pion momentum fluctuations imprinted at the phase boundary might be dissipated away, and thus 
"thermalized"to the final momentum space freeze-out, at about T = \OOMeV , as a pion experi- 
ences about 5 rescatterings on average during the final hadronic cascade [p^]. On the other hand a 
hadro-chemical K/norK/p ratio fluctuation would be preserved throughout the cascade as can be 
inferred from Fig. 

3.2 i^/TT Fluctuations 

We showed in Fig. ^ the first results concerning K/k (strangeness to entropy) fluctuations. 
We are looking for a reflection of the light and strange quark susceptibility maxima at /ig = 3 
shown in Fig. |3|, or of spinodal decomposition (Fig. The fluctuations are small from top SPS 
to RHIC energy in central collisions of Pb + Pb or Au +Au, but exhibit a remarkable rise toward 



the lowest SPS energies which can not be reproduced by UrQMD calculations [35]. The employed 
observable, Odyn, is based on a comparison of the widths a obtained for data and mixed background 
eventwise distributions respectively. We define 



■ sign{ol^,a - ^lix) ^\^La-^lix\ (3-4) 

the quantity that is illustrated in the right hand panel of Fig. ^. Note that in (4) crj^^^ and a^;^. are 
the relative fluctuation square width, i.e. ^data — (distributions )/< K/ti normalized to the 
ensemble data average of square eventwise K/n. Thus the % scale of Odyn- 
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Figure 16: Relative two particle pr correlator at SPS pS] and RHIC [620 energies. 



< K/n > represents the mean value of the data distribution shown in the left hand panel of Fig. ^. 
We note, for clarity, that this average ratio does not represent the global < K > /< tt > ra- 
tio, at given ^/s , which is obtained after systematic efficiency and acceptance corrections, from 
a completely different ensemble analysis. On the contrary, any event-by-event analysis is per- 
formed with the raw identified charged particle tracks, recorded in the individual event, without 
any corrections. Also, we repeat that the K/n ratio analyzed here [35] refers to both charge signs, 
+ /k^ + TT". Is the rise toward low ^/s (and thus toward /ig = SOOMeV) in the K/n fluc- 
tuation an indication of non-trivial dynamics [|^] at high /Xb? It has been shown [59, 6C] that Odyn 
in the K/n case where < K > is much smaller than < n > (note that these are the mean values of 
multiplicity of tracks in the acceptance of the analysis, not the corrected multiplicities in 4n or at 
—0.5 <y < +0.5) is given approximately by 



OdyniK/n) OC 



< {AK){An) > 
<K> 



(3.5) 



where the numerator expresses the true correlation effect but the denominator enters a feature 
characterizing the observational mode. This influence may dominate the signal as is shown in 
Fig. 18, the fall off in < ^ > thus causing the increase of the signal toward low ^ [^I]]. The 
same effect is exhibited by the STAR data at ^ = 62.4 and 200 GeV in Fig. |l9[ It shows 
the centrality dependence of a related fluctuation measure (see below), with a steep, monotonous 
increase toward peripheral collisions with decreasing accepted < K > and < n >. 
The fluctuation measure employed in the STAR analysis represents an alternative to Odyn employed 
by NA49: Vdyn expresses [53| the difference of the relative event multiplicities, — By 
definition the average value of this expression is zero. The corresponding variance of this quantity 
is 



K 



n 



<K > <n> 



<K^> 
<¥>2 



<Kn> 



+ ■ 



< 71^ > 



<i<:><7r> <n> 



2 ■ 



(3.6) 
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Figure 17: Transverse momentum eventwise fluctuation in central Pb + Pb coflisions [ |54| ] ; the SPS energies 
are represented by the corresponding values of jis- Also illustrated are enhancements expected from a critical 
point [||,p6l|7l Hi. 



For purely statistical fluctuations (Poisson statistics) this expression is reduced to 

^ + — ^ (3-7) 



<K> <n> 
and thus we have 

V -V V _ <K{K-\)> . <K7l> <7l{7l-\)> 



For further analysis of this variable see []63|]. It analyzes relative eventwise multiplicities like 
K/ < K >, similar to the approach taken [|^] with Odyn, and thus shares the feature to be inversely 
proportional to number of kaons (and pions) involved in the eventwise analysis, i.e. the < K >, 

< n > employed above. The STAR results in Fig. 19 for Vdyn {Kn), clearly exhibit this property. In 
these minimum bias runs the geometrical reconstruction volume (the acceptance) stays constant, of 
course, and unless there are systematic variations of the reconstruction efficiency with centrality, 
the accepted average eventwise multiplicities < ^ > or < tt > above are proportional also to the 
overall, "true" average midrapidity multiplicity < dN/drj > for charged particles, at each central- 
ity; this is assumed in employing dN/drj as the ordinate. We note that this procedure can not be 
extended to the NA49 data on energy dependence of K/n fluctuation (as is, in fact, done in ref. 
[|62|]) as the acceptance was changed with ^/s in the NA49 SPS runs, thus < K > and < tt > are not 
simply proportional to the variation of dN/drj with energy [54], in this case. These complications 
are, however, taken well into account in the fit shown in Fig. 18. 

We conclude that the rise in the K/n fluctuation signal represents, first of all, a consequence of 
considering the relative eventwise multiplicities and their fluctuations, both in a^yn and Vdyn- A 
new observable is required, to focus exclusively on the genuine correlation effect expressed by 

< {AK){An) >, or by the difference between < K/n > and < K > / < n >. The question of 
signals from the expected new physics thus remains open. 
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Figure 18: Relative fluctuation of + K /n^ + K , with (7^/,,„ in percent, similar as in Fig. 6. Also 
indicated is a dependence pQ] on the inverse number of average accepted kaons, at each 



3.3 Baryon - Strangeness Correlations: the K/p Fluctuation Signal 

A further fluctuation - correlation signal, addressing the connection between strangeness and baryon 
number, has recently been proposed by Koch, Majumder and Randrup []65[], as a diagnostic signal 
to study the QCD matter in the vicinity of the hadronization line in the [T, /x] diagram of Fig. |l[ The 
idea: in a medium consisting of deconfined quarks, the correlation between strangeness and baryon 
number is trivially caused by the strange quark quantum numbers. Baryon and strangeness number 
are strictly correlated. On the other hand, in a hadron-resonance state, strangeness is carried both 
by baryon number zero meons, and by hyperons carrying strangeness and baryon number. Here, 
the baryon-strangeness number correlation depends on the baryonchemical potential and, thus, on 
^/s of A + A collisions. 



A fascinating proposal arises: the baryon-strangeness correlation coefficient defined as [65] 



Cbs 



<BS> - <B><S> 
< 52 > - < 5 >2 



(3.9) 



is formulated to be unity in a simple quark-gluon "gas", with deviation from unity in a hadron- 
resonance system, depending on iig- On the one hand it is related to the quark susceptibihties, for 
u, d, s flavours ||6f 



Cbs ■ 



All 

zi 



ds 



(3.10) 



that have recently been obtained from a lattice calculation with three dynamical quark flavours at 
finite baryochemical potential. Please note that lattice QCD calculates susceptibilities as a function 
of (lattice) temperature, across the deconfinement temperature T^, such that in principle, Cbs lattice 
predictions are available both for the deconfined limit {T > Tc) and for QCD matter in the hadron- 
resonance limit (r < Tc). Thus, both the situations concerning baryon-strangeness correlation, 
anticipated in Cbs analysis, are thus accessible for the first time. 



Finally the BS correlation can be approximatively related [ |67| ] to the experimentally observable 
fluctuation of the eventwise kaon to proton ratio (more specifically: the to proton ratio). The 
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Figure 19: Centrality dependence of the relative K/ n fluctuation measure Vdyn in minimum bias collisions 
[|2|l of Am +Am at two RHIC energies. 



carries the major fraction of (anti-)strange quarks, the anti-hyperons not being accessible to 
any event-by-event analysis anyhow. The proton multiplicity at high pLs represents about one-half 
of the net baryon number density, the neutron remaining not identified. A preliminary theoretical 
analysis [ ]67| ] suggests an approximate relationship connecting the experimentally observed to 
proton fluctuation with Cbs'- 



< {AK+)^ > < {Apf > 4 Cbs 



thus closing the chain from lattice QCD at finite baryochemical potential, to observable K^/p 
fluctuation. 



We show in Fig. gO| the predictions made by Koch et al. [ |65| ] for an ideal parton gas, and for a 
hadron-resonance gas obtained from the statistical model along the chemical freeze-out line. Also 
shown are corresponding UrQMD calculations [^] which confirm the grand canonical predictions. 
Finally Fig. ^ shows Cbs from equ. (1 1) as extracted from recent 3 flavour lattice calculations [^], 
which indeed confirm at zero Hb the expected saturation at unity for T > Tc, and the decline at T 
well below Tc where hadronic freeze-out (of course not present in lattice calculations) would oc- 



cur. Such calculations are also becoming available [66] at finite Hb, such that the connection to 



/p fluctuation data could be explored. The first data refer to studies without kaon and proton 



charge separation, i.e. of {K^ + K )/p + p, and are shown in Fig. 21 along with the UrQMD 



simulation. The NA49 data [||] at SPS energies show a negative correlation except for the lowest 



energy where the signal shoots up. At RHIC energy the STAR data [ |68| ] give a positive correla- 
tion, and, with the notable exception of the point at 20 AGeV, the UrQMD results also exhibit this 
change of sign. It may be understood in terms of stopping net baryon number. At CERN energy 
< ^ > is much smaller than < p >, and the latter are thus not newly produced protons but the 
"valence" once shifted from beam to midrapidity. An upward fluctuation of the stopping mecha- 
nisms simultaneously leads to more protons at midrapidity and more stopped energy, available for 



meson production. In fact the n/p fluctuation measured by NA49 [ 35 ] at this energy also exhibits 
a negative correlation. At RHIC, on the contrary, both p and p at midrapidity are newly produced , 
competing with kaon pair production. 
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Figure 20: Baryon-strangeness correlation coefficient Cbs [p5|]. 

Left: expectation for a quark plasma, a grand canonical hadron-resonance gas, and for the transport model 
UrQMD 

Right: Cbs from equ. (11) extracted from a lattice calculation 1 27 1 at zero /Xb. 



The rise to positive fluctuation at the lowermost SPS energies remains unexplained. The fluctuation 
of /p, of key interest in the light of Cbs and lattice QCD analysis, will soon be available and 
may help to understand the low energy behaviour in Fig. 21. What we are looking for is a reflec- 
tion of the almost singular behaviour of lattice susceptibilities (Fig. ^ and ^ at /ig ^ 3 Tc, for light 
quarks but also for s quarks, to some extent. 



4. Outlook 

If it may be legitimate to write down a "belief" held by most if not all investigators we might 
express the expectation that the crossover nature of the phase boundary of Fig. |l| persists at least up 
to the domain of Hb = 200MeV, i.e. to top SPS energy. Our present interest concerns the physics 
in the vicinity of the phase boundary at higher Hb e.g. 

• a critical point of QCD, 

• an adjacent first order nature of phase transition, 

• a related " softest point" of the equation of state, governing eliptic flow, 

• a disentangling of the QCD chiral transition, from the confinement transition, 

• the related chiral restoration mechanism, 

• the exact position in the [T, Hb] plane of all this, 

• the existence, and domain, of a further QCD phase, and phase boundary, at high /ig ("quarky- 
onium"). 
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Figure 21: Preliminary results for the dynamical fluctuation of the + /T^ + ratio, from NA49 
and STAR [|8[. Also shown is the UrQMD model prediction. 



As a consequence of the expected need for substantial work at relatively low energies, say ^/s < 
20GeV, there arises an unfortunate split in the RHIC community, amusingly just coinciding with 
either accelerating or decelerating after RHIC injection from the AGS. It thus appears difficult to 
agree on an extended low energy scan at RHIC but this is exactly what is required. With the pos- 
sible exception of top SPS energy, ^/s =17.3 GeV, the statistics and systematics of the 2000-2002 
SPS low energy runs in the domain 250 < Hb < 450 MeV - which many consider the potentially 
most interesting - was far from adequate. Furthermore, the recent advances of lattice theory, both 
toward realistic quark masses and 3 flavour calculations, and toward higher statistical moments, 
and finite ^b, lead one to expect that all the topics enumerated above will come within reach of rig- 
orous theory. Of course, at this point we have to acknowledge, also, the present dichotomy among 
lattice QCD groups, about virtually everything - but then how about some experimental progress? 
For example an interesting physics topic, to be tackled next, consists in higher moments and cu- 
mulants of fluctuating quantities, such as quark densities, which exhibit interesting structures [27] 
in the vicinity of the critical (transition) temperature. These quantities are referring to the higher 
moments (beyond the width) of event by event fluctuating quantities, i.e. to symmetric (curtosis) 
or asymmetric (skewness) deviations from the normal distribution [27]. For example, Stephanov 
[ l69| ] has shown that such higher moments can better constrain the range ^ of the critical domain 
corresponding to a potential critical point (c.f. Fig. 17). First experimental data for the curtosis 
have been shown by STAR [70], for fluctuation of conserved charges at mid-rapidity. The impli- 
cation of conservation laws - the fact that a single A+A event in its totality represents a "grand 
microcanonical" ensemble - has been discussed in [jTll]. These considerations should give rise to a 
new class of physics observables. 
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